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^ ■ Abstract 

(— I I Lattice discretisation errors in the Landau gauge condition are examined. 

An improved gauge fixing algorithm in which 0{a^) errors are removed is 

|. ■ presented. O(a^) improvement of the gauge fixing condition improves com- 

\^ . parison with continuum Landau gauge in two ways: 1) through the elimina- 
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tion of 0(0^) errors and 2) through a secondary effect of reducing the size of 
higher-order errors. These results emphasise the importance of implementing 
an improved gauge fixing condition. 



I. INTRODUCTION 

Gauge fixing in lattice gauge theory simulations is crucial for many calculations. It 
is required for the study of gauge dependent quantities such as the gluon propagator [|l|, 
and is used to facilitate other techniques such as gauge-dependent fermion source smearing. 
However, the standard Landau gauge condition is the same as the continuum condition, 
J2fi d^A^ = 0, only to leading order in the lattice spacing a. 

There has been some study of alternate lattice definitions of the Landau gauge condition 
0. The focus of this letter is to use mean- field-improved perturbation theory 0] to compare 
different lattice definitions of the Landau gauge, and quantify the sizes of the discretisation 
errors. In particular, we derive a new O(a^) improved Landau-gauge- fixing functional which 
is central to estimating the discretisation errors made with the standard functionals. 

In section |I| we derive mean-field- improved expansions (in the lattice spacing and cou- 
pling constant) for three different definitions of the Landau gauge condition for the lattice. 
The use of an improved Landau gauge functional allows an estimate of the absolute error 
in standard lattice Landau gauge. The size of the error provides a strong argument for the 
use of an improved gauge fixing condition. 



II. LATTICE LANDAU GAUGE 



Gauge fixing on the lattice is achieved by maximising a functional whose extremum 
implies the gauge fixing condition. The usual Landau gauge fixing functional is 



where 



and 



= E ^Tr {U^{x) + u^ixy} , 



U^{x) = G{x)U^{x)G{x + 



G{x) =exp|-2E^"(^)^'' 



Taking the functional derivative of eqn (|2.1| ), we obtain 
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The gauge links are defined as 
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Ufj,{x) = Pexp |i(yf J dtAfj,(x + fit) 



(2.5) 



Connection with the continuum is made by Taylor-expanding Afj^{x+fit) about x, integrating 
term-by-term, and then expanding the exponential, typically to leading order in g, noting 
that errors are of 0{g'^a'^). Expanding eqn ( p.4| ), we obtain 



■G 



duj^'i^x) 



ga^ ^ Tr 



d^,A^{x) 
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a'dlA,{x) + —dlA,{x) + 0{a' 



rpa 



Oig'a 



(2.6) 



To lowest order in a, an extremum of eqn (|2.1| ) satisfies J2f_idfj_Afj_{x) = 0, which is the 
continuum Landau gauge condition. What it means in practice on the lattice is that 



j:[-j^dlA,ix)-n,j, 



(2.7) 



where Tii represents O(a^) and higher-order terms, as shown in eqn ( p.6|) . Naively one might 
hope that higher-order derivatives in the brackets of ( p. 61) are small, but it will be shown 
that the terms on the R.H.S. of eqn ( p.7| ) are large compared to the numerical accuracy 
possible in gauge fixing algorithms. 

An alternative gauge-fixing functional can be constructed using two-link terms, for ex- 
ample 



^2 = E ^Tr {u;:ix)u;:ix + a) + h.c.} 



(2.8) 



Taking the functional derivative yields 



bT^ 1 



(5a;'^(x) 



E Tr { [f/^^(x - 2ii)V';[x -li)- U^{x)U;:{x + fi) - h.c] T^} , (2.9) 



and expanding to 0{a^) we obtain 



6iu''{x) 



16 

d,A,{x) + y9jA,(x) + —a%A,{x) + 0{a^ 



+ 0{g^a 



(2.10) 



which again implies the continuum Landau-gauge-fixing condition to lowest order in a. 

(9(a^) errors can be removed from the gauge fixing condition by taking a linear combi- 
nation of the one-link and two-link functionals: 



'-^=^a^ETr{[9,A,(x)-— a^ajA,(x) + 0(c 



boj''{x) 



\ + 0{g^a^) 

(2.11) 
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where we have introduced the mean-field (tadpole) improvement parameter uq to ensure that 
our perturbative calculation is not spoiled by large renormalisations |Q. For the tadpole- 
improvement parameter we employ the plaquette measure 



Uo = {^TZeTr < Upi >y . (2.12) 



and note that the higher order g^a^ terms of eqns ( |2.6|) ( p.lO|) and (|2.11| ) are to be viewed 



in terms of the mean-field-improved perturbation theory [^. For future reference we shall 
define the improved functional, J^i^p = f-^f ~ 

Once a gauge-fixing functional has been defined, an algorithm must be chosen to perform 
the gauge fixing. We adopt a "steepest descents" approach [0. Collecting terms of 0{a'^) 
and higher into the Hi, we define 

Wo ^ 

= -2tga'j:l^d,A,{x) + ^dlA,{x)+n,^ , (2.13) 
= -2iga'Y.\d,A,{x) + '^&lA,{x) +'H^ , (2.14) 



Ainip(a;) = -Ai(a;) - -A2(a;) 

= -2^^a2ER^^(a;) +Hi^p}. (2.15) 

where the subscript, "traceless" denotes subtraction of the average of the colour-trace from 
each of the diagonal colour elements. The resulting gauge transformation is 

Gi{x) = exp |^Ai(x)| 

= l + -A,ix) + Oia'), (2.16) 

where a is a tuneable step-size parameter, and the index i is either 1, 2, or Imp. At each 
iteration Gi{x) is unitarised through an orthonormalisation procedure. The gauge fixing 
algorithm proceeds by calculating the relevant Aj in terms of the mean-field-improved links, 
and then applying the associated gauge transformation, eqn ( p.l6|) , to the gauge field. The 
algorithm is implemented in parallel, updating all links simultaneously, and is iterated until 
the Lattice Landau gauge condition is satisfied, to within some numerical accuracy. 
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III. DISCRETISATION ERRORS IN THE GAUGE FIXING CONDITION 



The approach to Landau gauge is usually measured by a quantity such as 

^^ = ^ETr{A,(x)A,(x)t} (3.1) 



which should tend to zero as the configuration becomes gauge fixed. Nc is the number of 
colours, i.e. 3. 

A configuration fixed using Ai(a;) will satisfy eqn (|2?7|). Substituting eqn ( |2.7| ) into 
eqn (|2.14| ) yields 



= -2zga' Y: l^^dlA,{x) - T^i + T^^j (3.2) 

and similarly, 

Aimp(x) = -2tga^ [^-^dlA^i^) - T^i + T^i^pj . (3.3) 

Since the improved measure has no O(a^) error of its own, eqn ( p.3| ) provides an estimate 
of the absolute size of these discretisation errors. 

IV. CALCULATIONS ON THE LATTICE 
A. The Gauge Action. 

The C(a^) tadpole- improved action is defined as 

Sg = ^E^eTr(l - U,,{x)) - ^5]7^eTr(l - f/,ect(a:)), (4.1) 

•J pi ^^^0 rect 

where the operators f/pi(x) and f/rect(a;) are defined 

Upiix) = U,{x)U,{x + fi)Ul{x + v)Ul{x), (4.2) 

and 

U,ect{x) = U^{x)Uy{x + jj)Uy{x + i) + jj)Ul{x + 2u)Ul{x + u)Ul{x) 

+ U^{x)U^{x + fi)U^{x + 2fi)Ul{x + fi + u)Ul{x + u)Ul{x). (4.3) 

The link product U^ccti.^) denotes the rectangular 1x2 and 2x1 plaquettes. Eqn ( [4.1| ) 
reproduces the continuum action as a — 0, provided that (3 takes the standard value of 
Q/g'^. 0{g'^a^) corrections to this action are estimated to be of the order of two to three 
percent 0. Note that our f3 = 6/g'^ differs from that used in Multiplication of our (3 

in eqn (ETI) by a factor of 5/3 reproduces their definition. 
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FIG. 1. The gauge fixing measures for a 6^* lattice with Wilson action at /3 = 6.0. This lattice 
was gauge fixed with Ai, so 6i drops steadily whilst 62 and ^imp plateau at much higher values. 



B. Numerical Simulations 



Gauge configurations are generated using the Cabbibo-Marinari pseudoheat-bath al- 
gorithm with three diagonal SU{2) subgroups. The mean link, uo, is averaged every 10 
sweeps and updated during thermalisation. 

For the exploration of gauge fixing errors we consider 6^ lattices at a variety of f3 for both 
standard Wilson and improved actions. For the standard Wilson action we consider /? = 5.7, 
6.0, and 6.2, corresponding to lattice spacings of 0.18, 0.10, and 0.07 fm respectively. For 
the improved action we consider /? = 3.92, 4.38, and 5.00, corresponding to lattice spacings 
of approximately 0.3, 0.2, and 0.1 fm respectively. 

The configurations were gauge fixed, using Conjugate Gradient Fourier Acceleration p| 
until 01 < IQ-^^. Oi^p and 62 were then measured, to see the size of the residual higher order 
terms. The evolution of the gauge fixing measures is shown for one of the lattices in fig. |1|. 
This procedure was then repeated, fixing with each of the other two functionals, and the 
results are shown in table |. The same procedure was performed with three (9(a^)-improved 
lattices, and the results are shown in table ||. 

Comparing eqn ( |2.7| ) with eqn (|3.3|), we can see that if we fix a configuration to Landau 
gauge by using the basic, one-link functional, the improved measure will consist entirely of 
the discretisation errors. Looking at table |, we see that at /? = 6.0, 6'ijnp = 0.058, i.e. 

^ETr|E^M^M(^) (e^^^^(^)) I = 0.058, (4.4) 
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e-, 


5.7 
6.0 

6.2 


0.865 
0.877 

0.88(3 


1 
1 

1 


0.0679 
0.0578 

0.0522 


0.611 
0.520 

0.169 


0.111 
0.111 

0.111 




f3 


Uo 


Functional 


^Imp 


01 


Simp 


5.7 

6.0 
6.2 


0.865 
0.877 
0.886 


2 
2 
2 


59.0 
65.1 
61.7 


33.2 
36.6 
34.7 


0.563 
0.563 
0.563 




P 




Functional 


Oi 


02 


Oi 

02 


5.7 
6.0 
6.2 


0.865 
0.877 
0.886 


Imp 
Imp 
Imp 


0.0427 
0.0367 
0.0332 


0.684 
0.588 
0.531 


0.0625 
0.0625 
0.0625 


TABLE I. Values of the gauge-fixing measures obtained usin 
lattices at three values of the lattice spacing, fixed to Landau g 
and improved functionals respectively. 


[g the Wilson gluon action on 6^ 
;auge with the one-link, two-link 


/? 


Uo 


Functional 




02 


^'imp 


3.92 
4.38 
5.00 


0.837 
0.880 
0.904 


1 
1 
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0.102 

0.0585 

0.0410 


0.921 
0.526 
0.369 


0.111 
0.111 
0.111 
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Uo 


Functional 


^Imp 


01 




3.92 
4.38 
5.00 


0.837 

0.880 
0.904 


2 
2 
2 


57.5 
53.4 
52.2 


32.3 
30.0 
29.4 


0.563 
0.563 
0.563 




P 


Uo 


Functional 


01 


02 


Oi 


3.92 
4.38 
5.00 


0.837 
0.880 
0.904 


Imp 
Imp 
Imp 


0.0638 
0.0366 
0.0261 


1.02 

0.586 

0.417 


0.0625 
0.0625 
0.0625 



TABLE II. Values of the gauge-fixing measures obtained using the improved gluon action on 
6^ lattices at three values of the lattice spacing, fixed to Landau gauge with the one-link, two-link 
and improved functionals respectively. 
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a substantial deviation from the continuum Landau gauge when compared with 9i < 10~^^. 

As a check of our simulations, we note that the definition for Ai^p, eqn (|2.15| ), provides 
a constraint on the measures when gauge fixed. For example 



'Imp _ I 12 

12 ^ 3 

Similarly, by fixing with A2(x) we expect 



6'imp { 12)^ = i ~ 111 (A 5) 



and fixing with Ai^p{x) leads to 

»i 1 



9, 16 (^-^l 



These ratios are reproduced by the data of tables | and 
A three-link functional can also be constructed, e.g. 



^sim = E ^Tr {U^{x - f.)U^{x)U^{x + /.) + h.c.} , (4.8) 



the functional derivative of which 



6uj-{x) 2 ^ 



leads to 
A3(x) 



E Tr { [t/i'(x - 3fx)U;:ix - 2fx)U^ix - /i) 

-^^ + + 2-") - h.c] T^} (4.9) 

-2z(7a2 E + l<^'df,A,{x) + ^a^M^) + ^ia')] + 0{g'a'). (4.10) 



If leading order errors dominate, then we should be able to make ratios like those above, but 
involving ^3. However, we find that the values taken from our simulations are very different 
from ratios based solely on leading, O(a^) errors. This indicates that the sum of higher-order 
errors is also significant. Whilst one could proceed to combine Ai, A2 and Ai^p to eliminate 
both 0(a^) and 0(0"^), it is likely that O{g'^o?) errors are of similar size to the 0{a^) errors. 
We will defer such an investigation to future work. 
A configuration fixed using Aimp(x) will satisfy 

E5^^^(a^)=E{-W- (4.11) 
Substituting this into eqn (|2.13|) yields 

Ai(x) = -2iga^ E | ^^lA.^x) +Hi- Hi^Jl . (4.12) 
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A comparison of this equation with eqn (|3.3|) reveals that the coefficients of the terms in 
curly brackets, expressing the discretisation errors in these two cases, differ only by an overall 
sign, which is lost in the calculation of the corresponding 9i. If the three different methods 
presented all fixed in exactly the same way, then the 9i^^ of a configuration fixed with Ai, 
would be equal to 9i when the configuration is fixed with Aimp- It is clear from the tables 
that they are not, signaling the higher-order derivative terms d^A^{x) take different values 
depending on the gauge fixing functional used. 

Examining the values in tables | and reveals that in every case 9i is smaller when 
we have fixed with the improved functional than ^i^p under the one-link functional. This 
suggests that the additional long range information used by the improved functional is 
producing a gauge fixed configuration with smaller, higher-order derivatives; a secondary 
effect of improvement. 

Equally, one can compare the value of ^2 when fixed using the one-link functional, and 9i 
when fixed using the two-link functional. In this case, their differences are rather large and 
are once again attributed to differences in the size of higher-order derivatives of the gauge 
field. The two-link functional is coarser, knows little about short range fluctuations, and 
fails to constrain higher-order derivatives. Similar conclusions are drawn from a comparision 
of ^2 fixed with the improved functional and 9i^^ fixed with the two link functional. 

It is also interesting to note that in terms of the absolute errors, the Wilson action at 
/3 = 6.0 is comparable to the improved lattice at /3 = 4.38, where the lattice spacing is three 
times larger. 

V. CONCLUSIONS 

We have fixed gluon field configurations to Landau gauge by three different functionals: 
one- link and a two-link functionals, both with 0{a?) errors, and an improved functional, 
with O(a^) errors. Using these functionals we have devised a method for estimating the 
discretisation errors involved. Our results indicate that order O{o?) improvement of the 
gauge fixing condition will improve comparison with the continuum Landau gauge in two 
ways: 1) through the elimination of 0{a?) errors and 2) through a secondary effect of 
reducing the size of higher-order errors. These conclusions are robust with respect to lattice 
spacing and we have also verified the stability of our conclusions by considering additional 
configurations to that presented here. We plan to investigate improved gauge fixing on larger 
volume lattices to see if these effects of improvement persist. Lattice Landau gauge, in its 
standard implementation, is substantially different from its continuum counterpart, despite 
fixing the Lattice gauge condition to one part in 10^^. 

ACKNOWLEDGEMENTS 

Thanks to Francis Vaughan and the Cooperative Research Centre for Advanced Compu- 
tational Systems (ACSys) for support in the development of parallel algorithms implemented 
in Connection Machine Fortran (CMF). DGR thanks Mike Peardon for helpful conversations 



9 



in the early stages of this work. This research is supported by the Austrahan Research Coun- 
ciL DGR acknowledges the support of PPARC through an Advanced Fellowship, and thanks 
FNAL for their hospitality during the course of this work. 



10 



REFERENCES 



[1] D.B. Leinweber, J-I. Skullerud, A.G. Williams & C. Parrinello, Phys. Rev. D 58, 031501 
(1998), [Eep-lat/9803015| ; ibid, submitted to Phys. Rev. D, |iep-lat/9811027| and references 
therein. 

[2] C.T.H. Davies et al, Phys. Rev. D 37, 1581 (1988). 

[3] L. Giusti, Nucl. Phys. B 498, 331 (1997); L. Giusti et. al, Phys. Lett. B 432, 196 
(1998). 

[4] G.P. Lepage & P.B. Mackenzie, Phys. Rev. D 48, 2250 (1993). 

[5] A. Cucchieri & T. Mendes, Phys. Rev. D 57, 3822 (1998). 

[6] N. Cabibbo & E. Marinari, Phys. Lett. B 119, 387 (1982). 

[7] M. Alford, W. Dimm & G.P. Lepage, Phys. Lett. B 361, 87 (1995). 

[8] F.X. Lee & D.B. Leinweber, Phys. Rev. D 59, 074504 (1999). 

[9] H.R. Fiebig & R.M. Woloshyn, Phys. Lett. B 385, 273 (1996). 



11 



